Band gaps of the primary metallic carbon nanotubes 
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Primary metallic, or small gap semiconducting nanotubes, are tubes with band gaps that arise 
solely from breaking the bond symmetry due to the curvature. We derive an analytic expression for 
these gaps by considering how a general symmetry breaking opens a gap in nanotubes with a well 
defined chiral wrapping vector. This approach provides a straightforward way to include all types 
of symmetry breaking effects, resulting in a simple unified gap equation as a function of chirality 
and deformations. 

PACS numbers: 71.2Q.Tx, 61.48.+C, 73.61. Wp 
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Recently, individual single wall nanotubes (SWNT) ex- 
hibiting small band gaps oL.the order of 10 meV where 
observed for the first timeE SWNT's can be classified 
according to their electronic band gap into three groups: 
semiconductors, small gap semiconductors and real met- 
als. A semiconducting gap arises when the graphite 
Fermi points are not allowed in the tube's Brillouin zone, 
which is given by distinct quantization lines according to 
the tube's circumferential boundary conditions. Such a 
gap is of the order of 1 eV and was predicted to scale with 
1/R where R is the tube's radiusJjQ a prediction that 
was later verified experimentally by a Scanning Tuniick 
ing Microscopy measurement of the density of states .Eftl 
The graphite Frami— points lie on a quantization line if 
mod( "~ m ) = 0,BtnI!J where n and m are the two inte- 
gers defining the tube's chiral vector. Tubes that sat- 
isfy this condition are called primary metallic. In a real 
nanotube, however, effects of curvature and deformation 
break the nearest neighbor bond symmetry, resulting in 
a shift of the two distinct Fermi points of graphite which 
lie at the corners of the hexagonal first Brillouin zone 
where the bonding and antibonding bands are degener- 
ate (If -points). This shift may open a gap depending 
on the position of the new if-points relative to the cir- 
cumferential quantization line. These gaps are about the 
value of room temperature and thus coined small gap 
semiconductors . 

Much interest has been devoted to the study of these 
small gap semiconducting tubes, which has led to a 
good basiCp-understanding. Already from numerical 
calculationsO it has been known for a while that only 
armchair tubes retain zero gap and therefore are truly 
metallic while chiral tubes open small gaps because of 
the intrinsic curvature. These gaps depend on the ra- 
dius and the chiral angle and were found to be lower 
than room temperature for tubes with radius exceeding 
6A. Other numerical studies also examined the effect. p£ 
structural deformation such as stressE£lE2l and twisrrJ'til 
on the electronic band gap of nanotubes where it was 
found that the armchair tubes open a gap under twist 
but are not affected by the uniaxial stress while an op- 



posite response arises injzigzag tubes. Analytic progress 
was made by Kane et alEj where the effects of deforma- 
tions were modeled as a perturbing vector potential in an 
effective mass Hamiltonian. Those calculations were able 
to verify the numerical findings for the gaps from the in- 
trinsic curvature as well as twists. After the completion 
of this work we also becapae aware of a more recent an- 
alytic work by Yang et alE3 where both twist and stress 
on the nanotubes were considered (but not the intrinsic 
curvature). _ . _ 

From the previous worksEm the gaps can m princi- 
ple be calculated for almost any shape of nanotube cither 
by using a numerical method or by determining the ex- 
act metric and curvature tensors. We now reconsider the 
effects of deformations on primary metallic nanotubes 
and formulate the problem in terms of a general sym- 
metry breaking in the tight binding model. This gives a 
straightforward analysis of the effect which results in a 
surprisingly simple and useful formula for the gap. Our 
compact expression combines the effects from intrinsic 
curvature, twist and stress and is only a function of the 
chiral wrapping vector (n, m) of the tube. This gives a 
quick and easy way to determine the gap and allows for 
a good insight in the physical effects as we will describe 
below. Our results are a direct consequence of any sym- 
metry breaking effect in the tight binding model of the 
graphite sheet in carbon nanotubes. 

Our starting point is the observation that the en- 
ergy separation between the bonding and antibonding 
bands according to the graphite tight binding scheme is 
27 53i=i e lk ' Ri , where Ri are the nearest neighbor bond 
vectors and 7 is the transfer integral which is the near- 
est neighbor Hamiltonian matrix element. E3 Since at zero 
temperature the bonding band is occupied and the an- 
tibonding empty, the Fermi points Kp lie at the band 
crossings, which are, for the unperturbed graphite, the 
six corners of the hexagonal first Brillouin zone. If we 
now break the symmetry of graphite and allow different 
transfer integrals 7, depending on the direction of the 
bonds Ri we arrive at a more general equation for the 
A:- vectors at which the bands crossEa 
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which defines the points of zero gap in fc-space K' F . 
For small changes % = 7 + Sji we expect small shifts 
in the band crossing location K' F = Kp + Ak. Since 
we are dealing with primary metallic tubes, we know 



that £ i=1 



AKfRt _ 



where Kp are the unperturbed 



Fermi points. Working in the nanotube coordinates (c, t) 
adopted from Ref. [l3] where c is the circumferential di- 
rection and i stands for the translational direction along 
the axis, the bond vectors are 
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where a 



2.49A is the length of the honeycomb unit 
vector and ct = Vn 2 + nm + m? is the circumference in 
units of a. Both inequivalent Fermi points in graphite 
give the same result when estimating the gap so it is suf- 
ficient to consider just one of them. For an unperturbed 
Fermi point we write 
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[(to + 2n)k c + mVSkt], 



(3) 



where k c and kt correspond to the fc-vectors along the 
circumferential and translational directions, respectively. 

In order to get the gap, we now need to know the dis- 
tance between the new Fermi point to the nearest quan- 
tization line at the quantized circumferential fc c -values. 
Since the quantization lines are parallel to kt, this dis- 
tance is given by Ak c the shift along the circumferential 
direction k c . Expanding Eq. ([!]) to linear order in the 
perturbations = 7* — 7 we find 

Ak c = j=[Sji(m — n) + <572(2n + m) — £73 (n + 2m)] 

(4) 

The gap is then obtained by exploiting the fact that close 
to the Fpfpii point the dispersion relation is linear and 
isotropicO 



Eg = V3aj\Ak c 



(5) 



We now want to determine the changes to the transfer 
integrals S^i due to the curvature and deformations. To 
first order this change is proportional to the change of the 
bond length between two neighboring carbon atoms, but 
may also be created by a misalignment of two neighboring 
7r orbitals. In general we find that we can always express 
the change of the nearest neighbor transfer integrals in 
terms of a bond deviation matrix D 



Sji = Ri • D • R % IR\ 



(0) 



This deviation matrix is useful for describing the effects of 
stress, twists, and curvature in a simple unified way as we 
will see below. All nanotubes have an intrinsic curvature 
which causes hybridization of the otherwise orthogonal 
7r and a orbitals. Since the a bands lie normally far 
from the Fermi energy, we only consider the hybridization 
effect on the it band, which crosses the Fermi point in 
primary matallic tubes as long as the tube's radius is 
R > 2.4A.y 

Following the calculations of Slater and Kostertj we 
can assume that the transfer integrals are proportional 
to cosine of the misalignment angle <f> between two neigh- 
boring 7r orbitals. (A calculation that takes into account 
the full rehybridization of all orbitals will be considered 
in a future study.) Using cos</> w 1 — | (^ _£ ) , we can 
immediately express the deformation tensor for the in- 
trinsic curvature in the basis of c and t 
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where a c _ c = a/\/3 is the carbon-carbon bond length 
and R = achjl-n is the tube's radius. From Eqs. (^-0) we 
find 
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This formula agrees remarkably well with previous nu- 
merical studiesEJ if we choose 7 = 2.5eV and also agrees 
with the results of Ref. [TJ|. Equation (JsJ) gives the en- 
ergy gap of all primary metallic nanotubes without any 
applied deformations. One observes that the armchair 
nanotubes to = n are the only real metallic tubes, while 
the primary metallic zigzag tubes (to = 0), open the 
highest gaps. 

Next we want to examine the effect of a general two 
dimensional structural deformation such as twists and 
stress on the gap in the primary metallic tubes. Our 
deformation can be written as R — > (I + S)^, where R is 
any vector on the tube's surface, I is the identity matrix 
and S is the deformation matrix 



S = 







(9) 



Here e c and tt are uniaxial stresses along the circumfer- 
ential and the translational directions and £ is the strain 
(nanotube twist). The bond deviation matrix is then 
given by D dcform = \R\bS, where b ~ 3.5eV/A is the 
linear change in the transfer integral with a change in 
the bond length 7^ — > 7, + 6|Ai2j|, and \R\ = a/y/3 is 
the bond length. We now use D dcform to obtain the 
of Eq. (||) and as we did with the curvature, inserting in 
Eq. (|) and using the dispersion relation we find 
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This equation is the response to a two-dimensional linear 
deformation within the graphite sheet. We notice from 
Eq. ([To]) that in the presence of equal uniaxial stresses in 
both directions c and t, a gap is not opened, as expected 
since the bonds would maintain their symmetry. The re- 
sponse of armchair and zigzag tubes is complimentary as 
noticed previously in the numerical studies, tia i.e. zigzag 
tubes have the maximum sensitivity for a uniaxial defor- 
mation (e c or €t) but are insensitive to twists £, while the 
opposite is true for armchairs, reaffirming that a twist 
deformation is the only possible source for an energy gap 
in the armchair tube. 

A uniaxial stress e c around the circumference corre- 
sponds to a global change of radius. This effect can 
come about as a time dependent deformation due to lat- 
tice vibrations in the breathing mode. Realistic static 
deformations on the other hand correspond to the intrin- 
sic curvature, a stress along the tube et and a twist £. 
Therefore the total gap from static deformations is given 
by combining Eq. (g) and Eq. @ with e c = 0. The total 
band equation now reads 
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which is the main result of our paper. 

In some cases it may be useful to express this formula 
in terms of the chiral angle a and the radius R of the 
tube, which gives 




FIG. 1. Energy gaps of primary metallic tubes of radius 
12 ± 0.5A with no twist £ = and a twist £ = 0.1% as a 
function of chirality. The continues lines correspond to the 
average radius (12 A) while the circles correspond to the ac- 
tual tubes. A very small twist can thus have a profound effect 
on the energy gaps. 



gaps have important consequences since they are gen- 
erally of the same order as room temperature for most 
primary metallic SWNT. Only armchair tubes are gener- 
ically gapless, but a very small twist induces a gap of 
the order of other small gap semiconducting nanotubes. 
Such a small twist, on the other hand, moves the gapless 
state to tubes with higher chirality. 
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In this fornuaur results are then consistent with previous 
calculationscJ where twists and the intrinsic curvature 
(but not stress) have been considered as a function of the 
chiral angle. After the completion of this work, a paper 
deriving the change in the band gap due to deformations 
was published Jl3 which is also consistent with the angu- 
lar dependence of the deformation part of Eq. ( |l2| ) (i.e. 
without the intrinsic curvature). 

We see that Eq. (|l^) shows an interesting interplay 
between curvature and deformation effects as plotted in 
in Fig. [I]. A very small twist can actually remove the 
gap due to the intrinsic curvature. For a given radius 
the gapless state is therefore moved towards tubes with 
higher chirality. 

In summary, we have presented a straightforward pro- 
cedure to calculate the energy gap induced by a general 
broken bond symmetry. This leads to a simple analytic 
expression for the band gaps from both the intrinsic cur- 
vature and applied deformations, which provides a quick 
and reliable way to estimate the physical effects. These 
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